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Table 1 

Summary of Lattice Results for OPE-Wg-improved continuum model. 



Interpolating Fields 


Mn 


\n 


e 


So 




(GeV) 


(GeV^) 




(GeV) 


Xi = £«'<= (m" Cj5 (f) m' 


0.938* 


0.016(3) 


10.0(3) 


1.97(12) 


X2 = e'''=(M''Crf')75M = 


Not seen 


0.0002(3)t 


2.43(12) 


1.84(9) 


(XiX2+X2Xi)/2 


Fixed 


0.0018(14) 


2.46(15) 


1.76(22) 


XsR = e'"=(«''C7^«')757''rf= 


0.96(4) 


0.035(6) 


6.6(2) 


2.18(12) 


XA = £"'<= («" Ct" m') «<= 


1.13(4) 


0.027(7) 


10.7(3) 


2.39(10) 


* Defines the lattice spacing a. 




t Inferred from 


XiXi and X1X2 


results. 



taminations in the correlation functions. How- 
ever, interplay between the fit parameters give 
rise to rather large uncertainties aX tj =9. 

4. ALTERNATE INTERPOLATORS 

These techniques are now used to investi- 
gate nucleon properties obtained from correlation 
functions of unconventional nucleon interpolating 
fields which become noisy prior to a clear ground 
state domination. The results are summarized in 
Table 1. 

There is a long history of argument over the op- 
timum nucleon interpolating field to be used for 
QCDSR analyses. The solution is now clear. This 
analysis indicates that the continuum model effec- 
tively removes excited state contaminations and 
allows the isolation of the ground state. Hence, 
one should choose interpolating fields such as XSR 
[6] that give good convergence of the OPE, as 
the corresponding increase in continuum contri- 
butions is tolerable. 

Ground state nucleon properties are indepen- 
dent of the interpolating field used to excite the 
baryon from the vacuum. This invariance is sat- 



Table 2 

Comparison with selected results for XsrXsr- 



Approach 



AAr(l GeV2) 
(GeV^) 



So 
(GeV) 



This work 


0.030(6) 


2.18(12) 


Ref. [7] 


0.032(1) 


1.92(5) 


Ref. [8] 


0.031(6) 


1.69(15) 


Ref. [2] 


0.022(4) 


<1.4 



isfied in a trivial manner. The interpolating field 
X2, which vanishes in the nonrelativistic limit, has 
negligible overlap with the nucleon ground state. 
The results of Table 2 for Ajv, evolved to a 
scale of 1 GeV^ in the leading log approxima- 
tion, compare favorably with other approaches. 
So is somewhat large as anticipated above. The 
smaller value of Ajv for Ref. [2] may be due to the 
omission of large a^ corrections to the Wilson co- 
efficient of the identity operator used to normalize 
Ajv- Such corrections could increase their result 
by 20%. 
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ratio of the continuum contributions as done in 
[2]. Our approach requires the use of tadpole- 
improvement [4] 



/2k 



1 



3k 

4:K„ 



1/2 



(10) 



to account for otherwise large renormalization 
factors. Remaining renormalization associated 
with composite operators is argued to be small 
[3]. The K dependence of these two wave function 
normalizations is very different, and is crucial to 
recovering the correct mass independence of the 
Wilson coefficient of the identity operator. 

With this approach, the effects of lattice 
anisotropy may be absorbed through a combina- 
tion of a larger continuum strength (^ > 1) and 
marginally larger continuum threshold (sq). To 
test the QCDSR method as closely as possible to 
its actual implementation, we fix the value for ^ 
obtained in a four parameter fit from t = 6 ^ 20. 
This value is similar to that obtained by fitting 
(7) to the first few time slices of the lattice data. 
The value of ^ depends on the choice of nucleon 
interpolating field and is largest for xn of (2) at 
10.0(3). This suggests the presence of significant 
anisotropy at the second time-slice following the 
source [3]. 

Infrared lattice artifacts are not a significant 
problem for this approach. The ultraviolet lattice 
cutoff is modeled in a manner similar to that for 
the continuum model. By the second time slice 
following the source, such corrections are found to 
be negligible. In the following we simply discard 
the source (t = 4) and first (t = 5) time slices 
from the fit of the correlation functions. 

3. NUCLEON CORRELATOR FITS 

Our purpose is to test whether the nucleon 
mass and coupling strength can be obtained ac- 
curately from a fit considering only the first few 
points of the correlation function. The lattice cor- 
relation functions [5] are fit with (9) in a three 
parameter search of Ajv, M^^ and sq, in analysis 
intervals from t = 6 ^ tf where tf ranges from 
9 through 23. The nucleon mass determined in 
each of the intervals is plotted as a function of tf 
in Figure 1. sq displays a similar plateau. 



3.0 

2.5 

S 2.0 
d 

II 

0.5 
0.0 



* i 



10 



12 



14 



16 
t. 



18 



20 22 



24 



Figure 1. The nucleon mass determined in each 
analysis interval plotted as a function of tf . In 
this and the following figure, bullets correspond 
to pole plus continuum fits from t = 6 ^ tf , and 
open squares illustrate a simple pole fit to the 
region tf —7 -^ tf , which is selected to give similar 
uncertainties in the nucleon mass at t = 20. 
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Figure 2. The nucleon coupling strength deter- 
mined in each analysis interval plotted as a func- 
tion oitf. Symbols are as in Figure 1. 



Figure 2 illustrates similar results for the cou- 
pling strength. It is interesting to see that the 
simple pole determination of Ajv fails to form 
a plateau at large time separations. Similar re- 
sults are seen for the larger values of k = 0.154 
and 0.156. The plateau in Mjv and Ajv for pole 
plus continuum fits indicates that the continuum 
model effectively accounts for excited state con- 
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interpolating fields. While there are formal field 
theoretic arguments indicating nucleon properties 
are independent of the interpolating field, it re- 
mains to demonstrate that this is in fact the case 
in practice. 

2. THE CONTINUUM MODEL 

Consider the following two-point function for 
the nucleon 



G(t,p) 



E 



■^tr 



r4(0|T{xivWxiv(0)}|0) (1) 



where 

XNix) = £"'<= {u'^''ix)Cj5d\x)) u'ix) 



(2) 



is the standard lattice interpolating field for the 
nucleon, and r4 = (1 -|- 74)/4 projects positive 
parity states for p = 0. At the phenomenological 
level, one inserts a complete set of states Ni and 
defines 



( I Xn{0) \ Ni,p,s)= \n,u{p,s) . 



(3) 



where the coupling strength, Ajv, , measures the 
ability of the interpolating field xn to annihilate 
the i'th nucleon excitation. For p = and Eu- 
clidean time t ^ CO, the ground state dominates 
and 



G(t) 



^2 - 



M^t 



The spectral representation is defined by 



Git) 



p{s) e "* ds , 



(4) 



(5) 



and the spectral density is p(s) = A|^ 6(s — Mn) + 
C(s) where ((s) provides the excited state contri- 
butions. 

While it may be tempting to fit the correla- 
tion function in the shorter time regime by in- 
cluding additional poles in the spectral density, 
such an approach fails in a number of ways [3]. 
The correlation function is probably best de- 
scribed by many states of diminishing coupling 
strengths and increasing widths. It may be that 
the QCDSR inspired continuum model is an effi- 
cient way of characterizing this physics. 



The form of the spectral density used in the 
continuum model is determined by the leading 
terms of the OPE surviving in the limit t ^ 0. In 
Euclidean space and coordinate gauge the quark 
propagator has the expansion 



^« " 27r2 x^ 



1 777, 



1 caa 



(27r)^ x'-" 



1 
"223 



QQ 



(6) 



and G{t) has the following OPE 



G(t) 



3-5V 1 28 m^a 14 mla^ 



287r4 



25 t^ 



25 t^ 



567r2 (: gg :) a^ 
"~75 t^ 



(7) 



The spectral density used in the continuum 
model is defined by equating (5) and (7). The 
continuum model is defined through the introduc- 
tion of a threshold which marks the effective onset 
of excited states in the spectral density. 



p(s) e "* ds 



(8) 



6 n— 1 ^ ^ 

EY^ — ^° r n 
Z^ L| fn-k ^"'-'" 



n = l k = 



where C„ and On are the Wilson coefficient and 
normal ordered operator of the term t~" in (7). 
The phenomenology of G(t) is summarized by 



G(t) 



\%e 



-Mjv* 



■c 



p(s) e "* ds . 



(9) 



Strictly speaking, ^ = 1 but here is optimized 
with Ajv, Mff, and sq to account for enhancement 
of the correlator in the short time regime due to 
lattice anisotropy. Ref. [2] found the anisotropy 
to be large for x — xq < 6 for free quark correlators 
and remain large in their interacting simulation 
at /3 = 5.7. At /3 = 5.9 there is some hope that 
anisotropy issues will be less problematic for the 
Fourier transformed correlators presented here. 
However, at very short times the quarks are es- 
sentially free and the anisotropy must be accom- 
modated. 

The nucleon coupling strength, Ajv, is deter- 
mined in absolute terms, without resorting to a 



OSU-PP #93-1121 



November, 1993 






in 

o 
o 

o 

-(— • 

I 

Oh 



A few points on point-to-point correlation functions 

Derek B. Leinweber^ 

^Department of Physics, The Ohio State University, 174 West 18th Avenue, Columbus, OH 43210-1106 

The short-time regime of QCD two-point correlation functions is examined through a QCD-Sum- Rule-inspired 
continuum model. QCD Sum Rule techniques are tested and alternate nucleon interpolating fields are discussed. 
The techniques presented here may be of practical use in determining heavy-light meson decay constants. 



1. INTRODUCTION 

In the quest for an ab tmtto determination of 
the low-lying hadron spectrum, lattice QCD in- 
vestigations have focused on the large Euclidean- 
time tails of three-momentum-projected two- 
point correlation functions. As such, the short 
time regime, where excited state contributions are 
significant, has simply been discarded. 

The lattice approach to QCD allows the de- 
termination of correlation functions deep in the 
nonperturbative regime. However, the QCD Sum 
Rule (QCDSR) method [1] is restricted to the 
near perturbative regime of the Operator Prod- 
uct Expansion (OPE). In this regime, one cannot 
ignore the contributions of excited states in QCD 
correlation functions. To account for these contri- 
butions, the so-called "crude continuum model" 
is introduced [1]. The contributions of this model 
relative to the ground state, whose properties one 
is really trying to determine, are not small. The 
validity of this model is relied upon to effectively 
remove the excited state contaminations. 

This investigation examines the physics in 
the near perturbative regime of point-to-point 
correlation functions where QCDSR analyses 
are performed. The QCD Continuum Model 
is constructed for three-momentum-projected 
Euclidean-time two-point functions following the 
techniques established in the QCDSR approach. 
This model will then be used as a probe of the 
physics represented in lattice QCD correlation 
functions, and as a test of QCDSR techniques. 

Some attention has recently been given to 
the behavior of lattice point-to-point correlation 
functions in the near perturbative regime. There 



[2] the focus is on space-like separated correla- 
tion functions. It was concluded that the QCDSR 
inspired continuum model was sufficient to de- 
scribe the lattice correlation functions over the 
calculated range. Using the entire lattice corre- 
lation function (including the deep nonperturba- 
tive regime), ground state masses were extracted 
and were found to agree with conventional lattice 
analyses. 

Here the emphasis is on determining ground 
state properties by examining only the first few 
points of the lattice correlation function. This 
is similar in spirit to QCDSR analyses. In this 
manner, the validity of the continuum model is 
rigorously tested. By extending the analysis in- 
terval of the correlation function deeper into the 
nonperturbative regime, the evolution of fitted 
ground state properties may be monitored. Con- 
ventional lattice results are recovered when the 
interval extends deep into the nonperturbative 
regime. A sensitivity to the analysis interval in 
the extracted parameters would indicate a failure 
of the QCDSR inspired continuum model. This 
approach also allows an examination of the im- 
portance of the higher order terms of the OPE 
in constructing continuum models. These terms 
were not investigated in [2]. 

One would also like to evaluate whether these 
techniques are useful in analyzing lattice QCD 
correlation functions. Such techniques may be 
of practical use for analyzing two-point correla- 
tion functions which become noisy prior to a clear 
ground state domination, such as heavy-light me- 
son correlators. 

Finally these techniques will be exploited to in- 
vestigate nucleon properties using unconventional 



